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A new method for computing the spin Hall conductivity for a two dimensional electron gas in 
the presence of the spin Orbit interaction is presented. The spin current is computed using the 
Many-Body wave function which is degenerated at zero momentum The degeneracy at K — gives 
rise to non-commuting Cartesian coordinates. The non-commuting Cartesian coordinate are a result 
of an effective Aharonov-Bohm vortex at K = 0. An explicit calculation for the Rashba model is 
presented. The conductivity is determined by the linear response theory which has two parts :A-a 
static spin Hall conductivity which is determined by the non-commuting coordinates and has the 
value jj^. B-a time dependent conductivity which renormalizes the static conductivity. The value 
of this renormalization depends on inelastic time scattering ,spin Orbit polarization energy and 
Zeeman energy. As a result the spin Hall conductivity vary between ^ and j^j-. 

In the absence of a Zeeman field we find that the long time behavior is given by the renormalized 
conductivity . 

For relative small magnetic field the Zeeman field allows to probe continuously the spin Hall 
conductivity from the static unrenormalized value ^ to the fully renormalized value . When 
the Zeeman energy exceeds the Fermi energy only one Fermi Dirac band is occupied and as a result 
the static Hall conductivity is half the static spin Hall conductivity. 

We compute the uniform magnetization without the Zeeman field and show that the spin current 
is covariantly conserved and satisfies effectively the continuity equation. 

The effect of a time reversal scattering potential due to a single impurity in the Rashba model 
causes the the spin Hall current to decrease with the size of the system. 



I. INTRODUCTION 

Two experimental groups 0,0 have attempted to confirm the spin Hall effect without reaching conclusive results. 
Experimentally one observes that spins of opposite sign accumulate on opposite sides of a semiconductor in response 
to an external electric field. The two experimental groups p| and have reached different conclusions. The Santa 
Barbara group 2j suggests that the spin accumulation is extrinsic and is caused by disorder; while the Hitachi group 
reports the experimental observation of an intrinsic spin Hall effect due to the spin Orbit interaction. Recently 
the Santa Barbara group Q have done a new experiment using a two -dimensional electron gas confined in the 
(110) direction in AlGaAs quantum wells with an electric field applied along the crystal axis. An out of plane spin 
polarization is reported which might be explained by the cubic Dresselhauss model ruling out the disorder extrinsic 
effect 0] . 

The proposed theories are also controversial, a side jump disorder induced model is proposed as a mechanism for 
the spin Hall effect 0, while other researchers have used the spin Orbit interaction to explain the possibility for a spin 
Hall current, [TL 1^. F7L I8L I9L IToL ITTL FT^ | ■ It seems that the definition of the spin current is also not clea r IT3 . The effect 
of static disorder introduces vertex corrections ^l| causing the spin Hall conductivity to vanish R, Illl fl^ . The 
conserved magnetic current in the presence of a orbital magnetic field has been investigated |10L Il5| and a proposal 
to verify the spin Hall effect by applying a magnetic field gradient and measuring a charge Hall current has been 
suggested 

Microscopically, the spin Orbit interaction emerges from the non-relativistic limit of the Dirac equation 0, ^} . At 
low energies the Dirac Hamiltonian for spin 1/2 electrons in an electromagnetic field is projected effectively into a two 
component Pauli Hamiltonian. This projection replaces the U (1) gauge fields — eA and —cAq by — eA — /ieo* x E 
and —eAo — • B where a is the Pauli matrix , \\ib is the magnetic moment, and -j-^b is the the Thomas 

precession [l^. E and B are the electric and magnetic fields. The Rashba Hamiltonian [ill [T2I IT7L l2t| used in Solid 
State physics is obtained after the replacement , ^^ ^jx E —* Hk so (ax £3) where £3 is the unit vector perpendicular 
to the two dimensional plane and E is the internal field which determines the strength of the spin Orbit momentum 
ksn = ^ 9 n 1 ) Ur\E\. For valence bands electrons the Rashba Hamiltonian is replaced by the Luttinger Hamiltonian 
GlIlllB. 
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The spin Hall conductivity has been computed in the literature 0, IToj without paying attention to the singularity 
at zero momentum. 

We construct an exact ground state in the presence of the spin orbit interaction which will be used to compute the 
spin Hall current. As an explicit example we investigate the Rashba Hamiltonian using periodic boundary conditions 
with the momentum restricted to a 2-torus T 2 . We use the spinor representation to compute the commutations 
relations for the coordinates in the momentum representation. This commutations emerge from the singular spinor 
transformation at K — used to diagonalized the Rashba model. The eigenstate spinor which solves the Rashba 
Hamiltonian is described by a singular vector field. By a singular point wc mean that the vector field is zero or is not 
smooth |2l| . In our case the spinor is multivalued at K = 0. The non-commutativity is a result of singular vector 
state at K = 0. In the geometrical language the non-commutativity is obtained from the knowledge of the connection 
|22|. The connection is computed in our case from the momentum derivative of the spinor-eigenstate \K ® ^ a (K) > 
,(a = 1,2 represents the two component spinor which is momentum dependent ) . The derivative of the connection 
allows us to compute the geometrical curvature which is interpreted as Aharonov Bohm effect in the momentum 
space with a fictitious non-zero magnetic field at K = 0. 

The theory which emerges from the vortex at K = gives rise to the non-commuting Cartesian coordinates . In 
order to compute the spin Hall conductivity we use the linear response theory. The conductivity for a constant electric 
field has two contributions: 

A-a, static spin Hall conductivity of determined by by zero momentum state. 

B- a time dependent conductivity which renormalizes the static conductivity which is a consequence of the gapless 
ground state . The time dependent corrections are given by the linear response term , -g^r[l — cos(2vps.k so t)] where 

vp s. = 2 ^ftj^ 30 ^ and fl = vp ,s.k so are the Fermi velocity and the the spin Orbit polarization frequency. 

The spin Hall conductivity is given by the sum of the two parts. We introduce the inelastic time scattering r s and 
compare it to the the Fermi Surface Spin — Orbit Polarization frequency f2 D = vp,s.k so . For times t — t s which 
obey fl a t < 1 we obtain that the spin Hall conductivity is given by the static vortex contribution For Q t s >> 1 

we perform a time average and obtain that the spin Hall conductivity takes the value, j^j- in agreement with the value 
reported in the literature 0, IS Ell EH E3 • According to ref. Q the spin Hall current is sensitive to the various 
scattering times . For time dependent electric fields our method allows for a simple separation between the static and 
the time dependent results. 

The Zeeman interaction breaks the SU(2) symmetry and modifies the spin Hall conductivity The presence of a 
weak Zeeman field allows the observation of a continuous variation of the spin Hall conductivity from the static value 
(obtained for Zeeman energy which are larger than the spin Orbit polarization energy) to the fully renormalized 

conductivity gjj- (obtained for a zero magnetic field). 

For extremely large Zeeman energy which are comparable to the the Fermi energy only one of the spin bands is 
populated. As a result we obtain that the static conductivity is j^J-. 

-The scattering effect of a non periodic time reversal invariant potential on the spin Hall current is investigated. 
We find that independent on the strength of the potential the spin Hall current vanishes with the size of the system. 
The scattering potential mixes the states. The zero momentum states remain degenerated in the presence of a time 
reversal invariant potential. The spectral function weight integral for the non zero momentum states diverges, as a 
result the eigenstate for the spin Hall current state has a vanishing overlap with the zero momentum state causing 
the current to vanish. Formally we find that the non-commuting Cartesian coordinates in the transformed basis (the 
basis which diagonalizes the scattering potential) obey the same transformation rule as the geometrical curvatures 
|2l|. The zero momentum states remain degenerated in the presence of a time reversal invariant potential. As a result 
the eigenstate of the scattering potential has a vanishing overlap with the zero momentum state. Therefore for a finite 
system the spin Hall current decreases with the size and vanishes for an infinite system. This result is not based on 
statistical averages and is in qualitative agreement with the transport results obtained in the literature 0, El 0, El • 

-The plan of this paper is as following. In chapter // we introduce the Rashba Hamiltonian using the spinor 
representation and compute the Many-Body wave function . In chapter III we use the spinor representation to 
represent the Cartesian coordinates . Using this spinor representation we compute the commutations relations. This 
commutations relations emerge from the singular spinor transformation which diagonalizes the Rashba Hamiltonian. 
In chapter IV we introduce the Heisenberg equation of motion needed for computing the spin Hall current. We 
define the spin Hall velocity using the Heisenberg and Interaction picture. In chapter V we compute the Spin-Hall 
current. The calculation contains two parts: A-The static spin Hall conductivity and B-The time dependent linear 
response spin Hall conductivity. In chapter VI we compute the uniform magnetization induced by the electric field. 
We find that for Q t s < 1 the uniform magnetization is finite and in the limit Q t s >> 1 it vanishes. In chapter 
VII we show that the spin current is covariantly conserved and the continuity equation is satisfied to first order in 
the electric field. In chapter VIII we add a Zeeman interaction to the Rashba Hamiltonian and compute the spin 
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Hall current. In chapter IX the scattering effect of a non periodic time reversal invariant potential on the spin Hall 
current is investigated. We have included three appendixes: In appendix A we replace the constant electric field by 
a time dependent vector potential. We derive the SU (2) Berry wave function using a time dependent momentum. 
In appendix B the SU(2) derivation of the spin current and the covariance conserved spin current are presented. In 
appendix C we compute the charge Hall current in the presence of both the spin orbit interaction and the magnetic 
field . 

II. THE MOMENTUM-SPINOR REPRESENTATION 

In this chapter we will compute the Many-Body ground state wave function for the Rashba model in the absence 
of an external field. We will show that this ground state has a vortex at K — 0. This wave function will be used 
to compute the spin Hall current. This method is applicable to a variety of spin orbit systems such as the cubic 
Dresselhaus model "'n"'type GaAs proposed to explain the spin Hall effect 0, 3- 

The Rashba model in the momentum eigenstate-spinor basis is given by, h = Iiq + h( ext ) where h$ is the Rashba 
Hamiltonian. 

~ h o = ^— [P- Tik so {a x e 3 )] 2 (1) 
2m 

The external potential is given by, ~h^ ext ) = —eE[ ext ^ ■ where e = — |e| is the electron charge, E[ ext ^ is the 

applied electric field at t > and R^ is the Cartesian coordinate in the % = 1 direction. A careful investigation 
of the two component Pauli Hamiltonian [T^ | suggest that the Rashba Hamiltonian in an external electric field has 
an additional external potential. This is seen when an external field is applied to the Pauli Hamiltonian which is a 
function of the total electric field, E = Ei int \ + E[ ext ^ . We substitute in the Pauli Hamiltonian the total electric 
field, E = % t) + E[ ext) and find ^-^ B a x (E (mt) + E[ ext) ) -> hk so (a x e 3 ) + (J^)^ x E[ ext) ) (E (lnt) is 
the absolute value of the internal electric field). As a result the external potential will contain an additional term, 
5h {ext) = -( m g fcs ° )o- ( ^p 2 E[ ext) . This potential is controlled by the magnetic charge and therefore in agreement with 

the literature will be ignored. The fluctuations of the spin Orbit are best described by the SU(2) Aharonov - Casher 
spin Orbit Hamiltonian |2^| (see appendix B). 

The Rashba Hamiltonian in eq.I is given in the momentum p = TiK and coordinate R — id/dK representation. We 
diagonalize this Hamiltonian using a spinor representation. The spinor which diagonalizes the Hamiltonian is singular 
at K = and therefore gives rise to non commuting coordinates. This methodology has been used for the Quantum 
Hall case where the magnetic Bloch functions have zero's in the magnetic Brillouin zone which gives rise to vortexes 
and to a quantized Hall effect [24[ (see eqs.2.f 3,3.7,3.8 and 3.9). For a solid with a periodic potential the plane 
wave functions is replaced by the Bloch functions u nt K{g)- The transformation from the plane wave representation 
to the Bloch representation modifies the coordinate representation from R = id/dK to f = id/dK + A n ^ n (K) where 
A n , m {K) = i J (27r) a u m K(o)~BT( u n,K{<T) H3- Zak has shown j2^ that for a degenerated point in the Brillouin zone the 
Berry phase A n ^ n (K) gives rise to a non zero band curvature, n) — -J^ x A n ^ n (K) giving rise to non-commuting 

coordinates, [A 1 " 1 , r' 2 )] = iQ(K, n). We will show that the eigenfunctions which diagonalizes the Rashba Hamiltonian 
have a singular point in the momentum space and therefore the methodology used in ref. t 25j applies to our case. 

The Rashba Hamiltonian has two eigenvectors, \K <g> ( a (K) >, a = 1,2 which replace the free particle-spinor 
\K <g> S >= \K > \S > where (|5 = 1 >) T = (1,0) , (\S = 2 >) T = (0, 1). Due to the periodic boundary conditions 
the vectors K^ and K^ are restricted to a two dimensional Torus, this ensures that for a finite system the state 
K = is included in the Many-Body wave function. The eigenvector spinor \K (g> ( a (K) > is obtained with the help 
of the SU(2) transformation used to diagonalize the Rashba Hamiltonian . We find, \K ® (, a {K) >—U((p(K),-d — 
f)\K ®S> with the SU{2) transformation which is given by, U{(p(K),tf = f ) = (e^-^f^^ \ (^ e ^n)^ {2) \ 

and tp(K) = arctan( ~^ T / ) is the azimuth angle in the plane, K = (K^ , K^). The Rashba Hamiltonian in the 
diagonal representation CT (3) fi_^M + £so has two eigenvalues, e a ( a ){K)) = o-(a) h fc ^° |g| + e so = 

~* h 2 k 2 ~~* 

e(\K\, er(a)),where a{a = 1) = 1 , o~(a = 2) = — 1 and e so = „ is the spin orbit single particle energy. At K = 

the eigenvalues are degenerated and the spinor \K C a {K) > is Singular (multivalued). The degeneracy at K = 
gives rise to a situation where the eigenvalue function which is a function of the momentum in the Brillouin zone 
returns to it's original value after a 47r rotation in the plane. Therefore the eigenvalue function has the topology of a 
connected sum of two tori T g= i . 
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Using the eigen-spinor basis we obtain the diagonal form of the Rashba hamiltonian, 

— J2z*(<x)(\K\)\K®UK) ><UK)®K\ (2) 

It is important to stress that the macroscopic spin current is determined by the ground state Many — Body wave 
function computed in the absence of the external field. 

For the Rashba model the wave function is determined by the products of the SU(2) transformations in the mo- 
mentum space which act on the state \G >= Y[g ® ^ > ' As a result we find that in the transformed frame 
the Pauli spin is aligned in the plane (the angle in the plane is determined by the value of the two component 
momenta). The Many-Body ground state \F.S. > for the Rashba model is determined by the product of the SU(2) 
rotations. 

\F.S. >= [] \K ® UK) >= H U(tp(K), $=~)\K®S> (3) 

K,a K,S 

The wave function \F.S. > is defined in terms of the two Fermi surface momenta ,K F and Kp. We introduce the 
notation , \F.S. >= \K F ,K F > to emphasize the fact that the occupation function is restricted by two maximal 
momenta K F and Kp. 



III. THE CARTESIAN COORDINATES IN THE SPINOR REPRESENTATION 



The basis \K®C^ a [K) > which builds the ground state wave function \F.S. > is used to find the Cartesian coordinates 
= 1,2 representation , 

r(i)= /(S/ 7^ii r( ^ p )\ K ^uK)><up)^p\ (4) 

with the matrix elements, r^\(K, P) — ® d ^ {i) 5(K — P) — \{a^) a ,i3 <8 S(K — P) d ^f^ , and momentum- spinor 
representation; 

= I^-j^m - ® = R (l) - ® where 1 is the identity operator and o-W is the Pauli matrix. 

The momentum derivative -j^jjj represents the Cartesian coordinate R^ for spinless particles and is the spinor 

connection which is a result of the SU(2) transformation which rotate the spin frame. 

The representation of the Cartesian coordinate in the momentum spinor representation given in eq.4 allows to 
compute the commutator for the Cartesian coordinates. In order to identify the singularity and to compute the two 
dimensional commutator it is advantageous to work in the complex plane, Z = + iK^ 2 \Z,= — iK^ 2 '. We 
perform a change of variables from the momentum plane K = (K^, K^) to the complex plane Z and Z. We find 
the commutator, 



[r«,r (2) ]^ (1) ^ (2) = — 

4 

The commutator in eq.5 is zero for K ^ 0. The right hand side of eq.5 represents the two dimensional delta function 
in the complex plane, 5^ 2 \K) — ~4=(i) = i^-( = ) which is a result of the multivalued phase tp(K — 0). At K = 

the eigenvalues are degenerated and the spinor \K ® Ca(K) > is Singular. The singularity in the plane at K = 
represent a vortex [2l| (see section 5.2 "'Closed Forms and Exact Forms ", page 156) and gives rise to non-commuting 
coordinates. The origin of this result follows from the fact that the energy spectrum forms a double Torus structure 
in the momentum space. The central result of this section is the result given by the commutator in eq.5. This result 
will be used to compute the spin Hall conductivity in the next chapters. 



^ dzdz -h^ )dzdz 



(5) 



IV. THE HEISENBERG EQUATION OF MOTION FOR NON-COMMUTING COORDINATES 

The singular transformation gives rise to the following commutation relations, [K^',K^] = 0, K^] = iSij 
and [rW, r^] = (1 — Si_j)^-a^2TTS 2 (K) replaces the commutator, [R^ , R^] = (6ij represents the Kronecker delta 
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function). Using the Heisenberg equation of motion we find that the time derivative of the momentum K^ 1 ' and 
the coordinate are linear functions of the electric field Ef xt ; ih dK> ' d ^ t ^ — —eEf xt and, V^^-fT, t) = ^—jr^- = 

(jf-){[r^ 2 \ ho] — eEf xt [r^ 2 ' , n 1 )]} is the velocity in the Heisenberg representation. Due to the non-commutativity of 
the coordinates ,we find that the velocity operator in the i = 2 direction depends on the electric field in the i = 1 
direction. The time dependent velocity Vi t H{K,t) (in the Heisenberg picture ) is given in terms of the Schroedinger 
Velocity Vi{K). The velocity and operator in the Schrodinger picture is given by, 

V itH (K,t = 0) = Vi(K) = ( dr( ' ) ]f ,f) )| f =o = (^)([r«,/To] - 6*rVV (1) ])|t=o (6) 

We will use the first quantized form to define the single particle operators in the momentum representation. For 
any single particle operator in the Schroedinger representation 0{K) we define the single particle operator in the 
Heisenberg Ojj(K,t) and interaction picture Oi(K,t). We will apply this formulation to the velocity and spin 
current. 

The time dependent velocity Vi^H(K,t) in the Heisenberg picture is given in terms of the Schroedinger 

Velocity V,{K) defined in eq.6. We find, V itH (K,t) = e^ h{R)t) V l {K)e^- fl(R)t = V it i{K,t) + 

(x) /o d *i[^,-r(-^'*)'- er j 1) (^'*i)] J5 i !Bt (*i) + ■■■ where Vij{K,t) = e* ft ~°^Vi(i?)eTrM^)* [ s the velocity in the 
interaction picture. 

The spin-velocity is defined according to the Noether's theorem [3(j (see appendix B). We will introduce the spin 
velocity in the Schroedinger picture , 

^\K) = \{V i {K),a^} + (7) 

{, }+ stands for the symmetric product and = U^(K)a^ A 'U(K) represents the transformed Pauli matrix with 
A = 1,2, 3. (This definition is equivalent to the symmetric projection of the 2x2 velocity matrix Vi(K) into the spin 
space and represents the velocity of the particle with a given spin polarization.) 

The spin-velocity in the i = 2 direction and polarization A — 3 takes in the Heisenberg representation the form, 

4%(K,t) = 4j(K,t) + (^) jf* dti^J^t), -erP&t^EFih) + ... (8) 

where ^(K^) and r { p{K,t) are defined according to the Interaction picture, 

4 3 )(iO) - e ^' ; °(^ t 4 3) (^)e*' l " o(/?)t = e!£f{K,t) +t ( *f xt - ) (K,t) (9) 

Where £^f\K,t) = es h "° W'^ 3,0 ) (j^)eTrM-fOt represents the homogeneous spin velocity in the Schroedinger 

picture /^'^(K) = \{jji{[ r ^ 2 \ ^o])|t=Oi &^}+- The second term in eq.9 represent the static spin-velocity generated 
by the non-commuting Cartesian. This term is generated by external electric field Ef xt at t = (the second term in 

eq.6). We find that the external spin velocity £^' ext '' (K) in the Interaction picture is the same as the Schroedinger 

picture £^ ext -\K) = £^f xL) (K,t) where ,tf ext \K) = f{^(-e£f* [r™ , r«])| t=0 , a^}+. 
The linear response result in our case is given by eq.8: 

The first term will give rise to the static linear response part , i£' ext ''(K) = i^f^ (K,t) — 

Mlk( — sEf xt [A 2 \ r^])\t=o, cr^ } + . This part is due to the non-commuting coordinates. 

The second term in eq.8 (the term which is linear in the external electric field Ef xt ) will give rise to the time 
dependent linear response when we substitute in the commutator of eq.8 the homogeneous spin velocity operator 

We will compute the expectation values using the second quantized form . We introduce the two component Spinors 
in the momentum representation ty>(K) and ^>(K) which act on the Many-Body ground state, \F.S. > . We introduce 
the Heisenberg operator in the second quantized form , Oh — J -^^s^(K)Os(K ,t)^ (K) . The second quantized 

(3) .... . 

form of the spin current H in the Heisenberg picture is obtained from the Heisenberg representation of the single 
particle operators. The spin current is obtained in the limit q — > by taking the expectation value of the spin current 
operator with respect the ground state \F.S. > . 
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4 3) = < F.S.\4%\F.S. > 

= Um^ {\) J ^<F.^.|vI/t(i?)(4 3 ) f (i?,t)e^(^))*(i?)|F.5.> 

ext. — static) ^ ^ time— dependent) (10) 

rn(K,t) is the coordinate operator in the Heisenberg picture and the exponential acts as a shift operator in the 
momentum space, e^ p( - R ^{K) = ^(K -q). 

The first term j( 3 - ext - statlc ) re p resen ts the static linear response current due to the the vortex at K = 0. The 
second term 5 j( 3 > tme dependent) . g ^ e n me dependent linear response . This term gives rise to renormalization effect 
of the static current. In the presence of a gap this renormalization effect are negligible. (This is the situation for the 
integer Quantum Hall where the conductivity is given by the static part.) 



V. THE SPIN-HALL CURRENT 

In order to compute the spin Hall current we will use the ground state wave function \F.S > character- 
ized by the Fermi — Dirac occupation function fp,D.[e(\K\,a(a)) — Ep] which at T = is given by the step 
function, 0[e(|i?|, <j(q.)) — Ep]. The spin Hall current has two parts ;the static linear response given by the first term 
in eq.8 and the time dependent linear response term given by the second term in eq.8. 

A- The static linear response current j( 3 ' ext -- statlc \ 

The formula for the static spin current in the i = 2 direction is given by the term, j( 3 > ext - statlc ) rp^g expectation 
value with respect the ground state \F.S. > at zero temperature gives, 

j^ext.-staUc) = ^ J ^^ <RS ^^jt)[£^\K)]^{K)\F.S.> 

= / /^^^P^^^^KiieuM)-^]^)^) 



(^ (1) ) 2 



d z { = )dZdZ - &z(^)dZdZ 

z z 
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= e -E^ t) - 7 ^Y\—l=-—l—\ = —E { r t) (11) 
2 1 2 2tt) 2 ^ l 2iri J Z 2wi J Z J 4tt 1 K ' 

a— 1 

The second row in eq.ll is a function of the commutator [r^ 1 ',^ 2 '}. The third row in eq.ll is obtained after we 
replace the commutator [A 1 ^,^ 2 ^] with his complex representation given in cq.5 . In the third row we recognize 
the two dimensional delta function ,8^ 2 \K) = i4=(-i) = \-§z~{=) which enable us to replace a two dimensional 
integral with a closed line integral in agreement with the Stokes theorem. The theorem allows us to replace the two 
dimensional integral over the full Fermi see with a closed line integral at the boundary of on the Fermi surface (the 
two Fermi surfaces, for spin up and spin down) . As a result we find that the static Spin — Hall current in eq.ll is 

replaced by a Fermi — Surface (the line integral in the fourth row) integral gives the exact value of ^ = ^ for 
the static Spin — Hall conductivity. The origin of this exact result, is caused by the current which is carried by the 
vortex at K = 0. The Fermi Dirac occupation function at K = takes the value of one for spin up and spin down, 
therefore the contribution from the two spin polarization in eq.ll are equal. 

B-The linear response time dependent linear response spin Hall current fij( 3 ' tlme ~ de P endent '> 

The linear response ,time dependent spin Hall current is given by the second term in eq.9. We use the interaction 
picture for the spin velocity , Cartesian coordinate and the Pauli matrix . £^f\K,t) — e^^ ^^ t £^2' \K)e~^ ^ s>t , 

Tr^f\K,t),(-e)r^(K,t)] = Tr[(^)K^ (K,t), ^(K^J^] 



-iO K(2) w^ Tr[ai w)'*/ wo] ( 12 ) 
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" 'Tr" ' represents the trace over the spin space . Using the interaction picture of the spin operators we evaluate the 
time dependent commutator, 

\af(K,t\af(KM)\ = (-2i)<r< 3 > sin[(^^)(t - t x )\ (13) 

m 

From equations 12and 13 we obtain the linear response time dependent term 5J^ ,tlme de P endent ^i w hich is computed 
from the second term in equation 8 by substituting i\ j\K, t) with l\ \ '(K,t). 

SJ (3,Ume- d e P en d en t) = J d?K ^ ^,0) ( ^ ^ ^(l) ( ^ fc)]^^) 

= f / d(^)[f^^+e so -^\K\-E F )-f F M^ + e so + ^\K\-E F )] 
57T J 2m 2m m 2m m 

■ [ d(-ti) sm[( 2MM )(t ~ t^E^ih) = ± f d(-2n o t x ) sin[(2fi (i - ti)\E?*(h) (14) 
Jo m 8?r 7 

At zero temperature the difference between the two Fermi - Dirac step functions for spin up and spin down is equal 
to the Fermi-Surface Spin-Orbit Polarization energy hfl Q multiplied by the delta function — h e so — Ep)]. The 
polarization frequency Vl — vp.s.k so is defined in term of the spin Orbit momentum k so = ^ 9 ~ 2 1 ' hb\E\ and the Fermi 
velocity v F .s. — 2 ^ et ^ ■ The momentum integration is replaced by the energy integration d( h 2 ^ ) restricted by 
the delta function <5[(%^ — h e so — Ep)] . As a result we find that the linear response time dependent current for a 
time dependent electric field is given by 8J^' ,tmie de P endent ^ = jl J* d(~2fl ti) sm[(2£l (t — ti)]Ef xt (ti) (see eq.14) 
For the special case where the electric field is time independent the time integration in eq. 14 gives the linear response 

Current, 5 J^Ume-dependent) = _e_fi _ C0S (2^ t)]Ef xt . 

Combining the two results ,the one obtained in equation 11 with the one given in equation 14 we obtain the spin 
Hall current for the time independent electric field Ef xt . 

j(3) _ j&ext.-static) + §J (3, time- dependent) = Zl E ^ + ± ^ _ c0s ( 2 Sl o t)}Ef xt (15) 

47T 8tt 

The result in equation 14 can be considered for different limits of inelastic scattering (the inelastic scattering time 
describes the temperature or other incoherent scattering such as spin waves or spin exchange) . The inelastic scattering 
is described by the time t s which has to be compared with the spin Orbit polarization frequency f2 Q . In the limit 
fl r s < 1 the second term in eq.14 can be ignored and one finds that the spin Hall conductivity is j^j-. In the opposite 
limit r2 T s >> 1 the second term can be replaced by a time average and we find that the spin Hall conductivity is 
given by jg- in agreement with the result obtained in the literature [ESS IH3 ■ ■ 

VI. THE MAGNETIZATION 

The application of an electric field in the i — 1 direction generates an uniform magnetization with an in plane 
spin polarization M^ A=2 ^ (in the i = 2) direction which is time dependent and vanishes in the long time limit. The 
uniform magnetization is given by , 

M (A) {t) = *Jj0_ < RS .\yt(K)a^\K,t)*(K)\F.S. > (16) 

As a result of the SU{2) transformation the spin operators are transformed; a^^K) = cosip(K) — sinip(K), 
<j( 2 ) (K) — c' 2 ' cos ip(K) — ct*- 3 -* smip(K) and — — cr' 1 '. The Pauli matrices A = 1, 2, 3 in the Heisenberg represen- 
tation are represented in terms of the Pauli matrices and the coordinate r^p in the Interaction picture. 

c4 A) (i?, f) - o-\ A) (K, t)-~J (K, t), -er ( p (K, t^E^ih)... (17) 

Using the explicit form of the Pauli matrices and the SU(2) angular <p(K) dependence, smtp(K) = — -j^-pand 
cos <p(K) = i y^r we find that the only non zero term is M^ 2 \ This is the magnetization with the A = 2 polarization 
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(this is the only term which is invariant under the transformation K — * —K and has a finite trace over the Pauli 
matrices ). We have, 

M^(t) = ~J 0y<F.S\&{K){ C osv{K)[af\K,t)-~j*tt 

- sin^(#)[<xf (K,t) - 1 jT dt 1 [af\K 1 t)-erf\KM)]ET t (t 1 m(K)\F.S. > (18) 

The commutator in the last equation is proportional to d J^i] — ~ j^jj therefore only the product with the term 
proportional to cos tp(K) = will be nonzero. As a result we obtain, 

7L,-iw-e, /" d 2 iT 



M(2) W ^ <5><T><T> / 



(2tt)2 



— e 



,ft 2 K 2 h 2 k SOt ^ t , ^ 2 K 2 h 2 k SOt ^ t 



Jo "I 



' dK( —)S(— +e so -E F ) dt lC os( — (t - hfiE^h) (19) 



8tt J m 2m 

For a constant electric field we perform the time integral and find the time dependent magnetization. 

where S1 Q is the the polarization frequency fl = VF.s.kso defined in the previous chapter. We observe that for inelastic 

scattering times r s which obey f2 T s < 1 the magnetization is finite, w (|)(y^)fc so [ ~ e I — — ]. In the opposite 

limit we perform the time average and find that the uniform magnetization vanishes. 

VII. THE COVARIANTLY CONSERVED SPIN CURRENT 



Following Noether's theorem |3jJ,|31| (see Appendix B) we compute the spin currents . We apply this formulation 
to the Rashba model and find from cq._B9 that the spin current is covariantly conserved. The three spin polarizations 
currents obey in the long wave limit limit q — » the following continuity equation ; dtJ$ + iq ■ J^(q) = 0, 
dtJ(p + iq - J^ 2 \q) = 0. Only the A = 3 component violates the continuity equation 

d t 4 3) + iq- J^(q) = 2k so (j[ 1) (q) + J (2) (q)) 

^ 2k so {\)<F.S.\ J |^^(J?,f)[i{^(^ 

,7k f d 2 K „„ lT _ L ,-, r l, ( 3 ),-, ,,-t .s K 



= 2fc so (-) J — <F.S.\^+(K)[-{a^(K,t),(V H (K,t) x -^)}+ 

+ ^ 2) (K,t), • 4-)}+]*(^)|F.5. > (21) 

This equation has been obtained from the Heisenberg representation of the spin velocity and the SU(2) transformed 
Pauli matrices. ( — cosip(K) — smip(K),a^ — cosip(K)) — sinip(K)). The right hand side of 
equation 21 is evaluated using the Heisenberg representation of the operators. The velocity and spin operators are 
a function of the external electric field, see eqs.8and 17. To first order in the electric field we observed that the 
right hand side term vanishes. (The first term on the right hand side of the equation vanishes as a result of to the 
momentum integration .and the second term vanishes as a result of the expectation value with respect the Pauli 
matrices.) Therefore we conclude that the continuity equation is effectively satisfied. 
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VIII. THE SPIN HALL EFFECT IN THE PRESENCE OF A ZEEMAN INTERACTION 

From the analysis in chapters IV and V we learn that the static spin current is determined by the state K = . 
Since the Fcrmi-Dirac step function for the state K = at zero temperature is not affected by a Zecman field (if this 
energy is much less than the Fermi energy) the static spin-Hall current will be same as for the zero magnetic field 
(see eq.ll). The effect of the magnetic field will be to renormalize the static current. 

We will present first the modification caused by the Zeeman field to the Rashba model. The Zeeman energy 
&3 = ^gusB ( B is the effective magnetic field) changes the the single particle eigenvalues from , e„( a ){K)) = 

n*gf a{a) ^Ml + tso t0 €a{a){K)) = _ a{a)A (K) + e so where A(K) = y/ % + (^^. The SU(2) 

transformation is replaced by , U{Lp{K),-d{K)) = e (-V2)y(i?)<x (3) e M/2)#(i?)o-< 2 > with the polar angle = | ( zer0 

magnetic field) replaced by tan $ (if) = h (finite magnetic field). 

Following the steps described in equations 3 — 6 we find that the polar angle d(K) modifies the transformed the 
Cartesian coordinates i — 1,2. 

rc«(f) . - y> ^m> d H§ - - ™m> d i§) w 

and the commutation relations are transformed , 



r (l) r (2) 



cr (1) sin?9(if) + cr w cos??(if) 2n5' 2 (K) 



r(3) 



(23) 



where the transformed Pauli matrix is given by, — > Wa^U — cos -d{K) — sin d(K) and the delta 
function ,6W(K) = \^{\) = ±&(±) 

At zero temperature the ground state is described by the two Fermi — Dirac for spin up and spin down)step 
functions, 8[e(\K\, a(a)) — Ef\- We will consider two cases: 

a) The magnetic Zeeman energy is much less than the Fermi energy therefore the sate K — with spin up or spin 
down have equal occupation, 8[e so + b 3 — Ep] = 9{—b 3 + e so — Ep] = 1- For this case the static spin Hall current is 
identical to the case without the Zeeman field .Following the method used in section V we find that the spin Hall 
current is given by two contributions ;the static one (see eq.ll) and the the time dependent contribution (see eq.14). 
The sum of the two parts is given by an expression similar to the one given in eq.15. 



4" - + ^(^0^J)[1 - ™(2«t + (24) 

The static spin Hall conductivity is determined by the degenerate state at K — and is given by the first term in 
eq.24 which is identical with the first term in cq.15. The second term in eq.24 represents the time dependent linear 
response part. Here we observe that the Zeeman magnetic interaction controls the renormalization effects. The Zeeman 

interaction modifies the value of the commutator given in cq.12. which is replaced by the term sin 2 d(K) = ( (nny+b^ )■ 
When the Zeeman energies b% is larger than polarization energy TiSl the linear response term can be ignored and the 
spin Hall current is determined by the static solution, ~ j^E^ xt \ For Zeeman energies which are comparable 
or less than the polari 

is given by, M[l - §(faj o)2+ &g n 

Therefore in this case the magnetic field can be used to vary continuously the spin Hall conductivity from the 
static conductivity ^ for A- >> 1 ( for large magnetic fields comparable to the polarization energy ) to the fully 

renormalized spin Hall conductivity j^j- for 6 3 = (zero magnetic). 

b)-When the Zeeman interaction is comparable or larger than the Fermi energy only one of the Dirac step function 
contributes. For this case we find a spin Hall conductivity is determined by the static part (without normalization 
by the time dependent linear response term). In this case only one of the Dirac step functions contribute to the 
conductivity 



or less than the polarization energy we find in the long time limit, (fi 2 , + (^f) 2 ) 2 ^ >> 1 that the spin Hall current 

|e| r , _ if (hn f xi 
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IX. THE SPIN-HALL EFEECT IN THE PRESENCE OF A NON-PERIODIC SCATTERING 

POTENTIAL 

We will investigate the effect of a non-periodic time reversal invariant potential V(r) on the Rashba hamiltonian 
in the absence of inelastic and scattering and Zeeman field. We will see that the spin Hall current depends on the 
ratio where L is the size of the system and A so is the spin orbit length defined by the inverse of the spin Orbit 

momentum A so = |^ where k so = = 2=L/j, B \E\. 

The hamiltonian, h = ho + V(r) is represented with the help of the the Rashba eigenstates basis, 

h = f 0^^e a{a) (\K\)\K®UK)><UK)®K\ 

+ I JTX2 I E E P ^ K >< ^ p l ( 25 ) 

J V / J Q=1/J=1 

Where the matrix elements of the scattering potential are given by the product of the Fourier component in the 
momentum representation, V(7?-P) with the SU(2) matrix elements, U{Lp(K),-d = f ) = (e^/ 2 ^^ 3 ^ ( e ( _i / 2 )i <72 ) . 

We define the transformed potential, V a ,p{K - P) = Ea=i U*(<p(K),0 = %) a ,\V(K - P)U((p(P),0 = f )a,/3- The 
eigenstates \$( a \K) > of the hamiltonian, h = h + V(r) replaces the eigenstate \K ® Ca(K) > given in eq.2. 

1 2 

'( K ) A=l 

+ E EA Q) (^A)|P®CA(n>] (26) 

The new eigenstate \&( a \K) > is given in terms of the amplitude A^(P, A) which represents the overlap between 
the eigenstate \$ a (K) > and the spinor, \P (g> (,\(P) >■ The amplitude \<f>( a \K) >is computed with the help of 
the matrix elements of the scattering potential given in cq.26 .The function N^ a \K) represents the normalization 
factor defined by the orthonormality condition, < $( a \K)\&P\P) >= 5 a jj5gP. The non-periodic potential is time 
reversal invariant therefore the zero momentum eigenstate \& a \K = 0) > remains double degenerated! In order 
to deal with the state K = we replace the integrals over the momenta with a discrete sum and replace the two 
dimensional delta function 5 2 {K = 0) with the Kronecker delta function 6g Q . We will use the new basis given in eq. 
26 to represent the coordinates (similar to eq.4). We find the matrix elements of the Cartesian commutator in the 
new basis (the basis due to the scattering potential). 

< ^ a \K)\[r^,r^]\^\P) >= Hl(<y {1) Up}S R p S R o [N^\0)N^(0)}-i (27) 

Using the result given in eq.27 we compute the spin Hall current following the steps given in cq.ll. The static current 
is determined by the normalization function N^(K = 0) given by eq.26. 

2 2 
T (3,ext. static) _ S „(ext) 1 V"^ A r 1 1 _ r 1 \ 1 ~ e p( ext ) (OS\ 

a — l K a—1 V / 

The normalization factor N^(K = 0) is a function of the amplitude A^L Q (P, A) (the projection of the eigenstate 
= 0) > on the spinor \P eg) Ca(-P) > ). From eq.26 we determine he normalization factor, 

The amplitude A^L {P, A) is a function of the matrix elements V(q) and is obtained within perturbation theory , 
\A^JP, a)\ 2 = \ , K -i V ^ K ~ Q ~ pS> 1 2 . This allows to compute the current in eq.28. We substitute the amplitude 
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A^^ Q (P, A) into the normalization function N^ a \K = 0) and obtain from eq.28 the current; 

(3,^.-^0 = [1+ rn/h 2 1 rl[ ^ t)] (30) 

2 1 (h 2 /2m)k 2 so 2tt g(l + i g 2 )2 47r v 7 

In eq.30 we have introduced the dimensionlcss momentum q (the momentum has been normalized by the spin orbit 
momentum k so ). The last integral is evaluated with the help of the dimensionless infrared cutoff K ™ in = where 
L is the size of the two dimensional system, A so is the spin orbit wavelength and A is an arbitrary uftraviolet cutoff 
(which does not affect the integration). The momentum integration can be performed when the matrix elements 
(V(q)) 2 are replaced by a momentum independent potential. This introduces the scattering strength exponent 7 SC 

= (V(q)) 2 m/h 2 
' < sc — h 2 /mk 2 a ■ 

We observe that for an infinite systems the integral diverges causing the spin Hall current to decreases with the 
size of the system L, 

j^t.-staUc) = [1+7scLn( ^ r l [ M S (-*)] 

Next we add the time dependent part spin Hall current. This is done using equations 10 — 15. In the absence of 
inelastic scattering f2 T s — » oo we find following the result in eq.15 that the the total spin Hall current is given by, 

j^li + vMA)]-^^), (31) 

A so tin 

We observe that for an infinite systems the integral diverges causing the spin Hall current to decreases with the 
size of the system L. The spin Hall current vanishes for an infinite system since the scattering causes the weight 
of the zero momentum state to spread over to all other states. Since the integrated weight for the non zero states 
diverges in two dimensions we find that the projection of the eigenstate with the K — state vanishes. This result has 
been obtained using general arguments, the only assumption is that the absolute value of the square of the scattering 
potential is momentum independent, such a situation is realized for a single impurity potential. For many impurities 
we obtain the same result as for the single impurity if we replace the square of the scattering potential by an ensemble 
average over the impurity configuration, V(q)) 2 «< (V(q)) 2 > configuration-average— V 2 C and obtain the scattering 

exponent 7 SC = jfij^j^r- For many impurities the fluctuations of the scattering potential will cause the current to 
vanishes faster than our power law prediction. 



X. CONCLUSION 



A new formulation for the spin Hall current based on exact ground state of the Rashba Hamiltonian has been 
presented. This formulation shows that the degeneracy of the wave function at K = gives rise to non-commuting 
coordinates. Due to the fact that the spin Orbit wave function has no gap the spin Hall conductivity is composed 
from two contribution ,a static part caused by the non-commuting coordinates and a time dependent linear response 
part. When the Zeeman interaction is larger than the polarization energy the spin Hall conductivity is determined 
by the static part. The uniform magnetization is computed and we show that the spin current which is covariantly 
conserved is effectively conserved to first order in the electric field. 

Using this new formulation we find the exact value for the spin Hall conductivity which vary from to 
An exact derivation for the vanishing the spin Hall current caused by a scattering potentials has been presented. 
The spin Hall current vanishes for an infinite system since the scattering causes the weight of the zero momentum 
state to spread over to all other states. Since the integrated weight for the non zero states diverges in two dimensions 
we find that the projection of the eigenstate with the K = state is zero. This result has been obtained using general 
arguments, the only assumption is that the absolute value of the square of the scattering potential is momentum 
independent, such a situation is realized for a single impurity potential. 



APPENDIX A: THE SU(2) BERRY PHASE 

The effect of the static electric field is equivalent to a time dependent vector potential . As a result obtain a time 
dependent momentum vector, K^\t) = K^(0) + j-E[ ext h, K^(t) = K^(0) .Therefore according to the adiabatic 
approximation we replace the spinors \K ® C, a {K) > by the time dependent state, \K(t) ® ( a (K(t)) >. According to 
Berry [HISl one nnds instead of the spinor \K(t) <g> £ a (K (i)) >= cxp[^- f Q e a{a) (K{t')) dt']\K(0) ® ( a (K(0)) > the 
SU(2) time dependent solution \^ a (K, t) > is given by, 
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t • ,* 



V a (K,t)>=Texp[-i n(K(t)dt]exp[— e <a) (K(t))dt]\K(0)®t a (K(0))> (Al) 



o 

We observe that the Berry phase has been replaced by a time ordered matrix, Texp[— i J Q * Q(K(t ) ] with the 
SU{2) matrix elements, to a , a (K(t)) =< K(t) ® C a (if(*))| ^^(i) <S> ( a (K(t)) > and non-diagonal matrix elements, 

n a=1 , /3=2 (^(t)) = exp[^ f '(e CT(Q=1) (tf(t'))-e ff( £ =2) (tf(t^ >• Thc 

SU(2) spinor l^^if, t) > is used to construct the time dependent representation for the coordinate and the velocity 
operator. 

2 2 



T(e l fo n( - K (*' ))dt ') e % /o £Q(if(t ' ))rft 'r W /3 ( J ft:,P)e^ i /o £ " (K(t ' }) dt ' T{e~ l /o' !2(if(f ' ) ' lt ') (A2) 



where r^o(Jf, P) is the matrix element given in eq.4. The time derivative of the coordinate representation given in 
eq.Al represents the velocity operator. 

dr(2) <™ } = l[r( 2 )(^(t)),/ l0 (^))] + i[r( 2 )(ii:W),OTO))] 
dr in in 

= ±;[rW(K(t)),h (K(t))] + U-e)E^ [r™ (K(t)), r« (K(t))] (A3) 
in in 

The first term in eq.^42 represents the velocity in the absence of the external field. The second term describes the effect 
of the external field proportional to the commutator of the Cartesian coordinates. This can be seen in the following 
way, the matrix Q(K(t)) is equal to the product between the time derivative of the momentum and a new matrix 
( dirt 1 ) )• For this matrix we replace the momentum derivative by the time derivative momentum derivative. We 
find, Q(K(t)) — dK ^ t - Cl( -jSjj) )■ The momentum time derivative is proportional to the electric field and the matrix 

■rSfT) is identified with the Cartesian coordinate ,[r^ 2 \ Vl{ d d {1) )] = [A 2 \A^]. As a result of the SU(2) adiabatic 
evolution , the velocity operator given in eq.A2 is equivalent to the the Heisenberg equation of motion given in eq.6 
(with the same commutation rules as given in eq.4). 



APPENDIX B: GAUGE INVARIANCE IN THE MOMENTUM SPACE-THE COVARIANTLY 

CONSERVED SPIN CURRENT 

In this appendix we will the consider the U(l) x SU(2)( spin } gauge invariance of Ahharon-Casher (A-C)[2j| model 
for a periodic lattice in the momentum space. The A-C effect is viewed as a moving magnetic moment which is 
equivalent to a moving magnetic current, due to the Lorentz transformation an induced electric-charge appears in the 
laboratory frame, which interacts with an electrostatic potential. 

H = Ye A{r) ~ {g ~2hc^ {& X ^ (r))]2 " I - { ^ XI) (B1) 

where A(r) is the U(l) electromagnetic vector potential, A^(r) = (a x E(r)) is the SU(2) "electrostatic vector 
potential" (a is the Pauli matrix and E(r) is the electrostatic field), /is is the Bohr magneton, g ss 2 is the gyro 
magnetic factor with g — > (g — l)/2 being the Thomas precession and "e" is the electrostatic charge. 

For a periodic lattice we use the quasi-momentum representation with momentum integration restricted to the 
first Brillouin zone. The single particle energy and the two component spinor obey the symmetry in the momentum 
space (for simplicity we replace the Bloch wave function by the free particle representation), e(K) = e(K + G) , 
^> + (K) = ^ + (K + G), ^(K) = ^(K + G) where G is the reciprocal Lattice vector. In order to describe the two 
dimensional Rashba model we take the z component of the electric field to be constant in space . We introduce the 
static SU(2) Rashba term, 
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J2a=i,2^a(K)<7^ = k so (a x e 3 ) 

Only two components are nonzero; ,a=2) {K)<7^ a=2 ^ = — k so a^ 2 \ lo(ja=i){K)o-^ a=1 ^ = k so a^ and 
uj(i t A=i)(K)a^ A=1 ^ = ^J(2,a=2)(K)<7^ a=2 ^ — 0. The fluctuation of the electric field gives rise to the SU(2) gauge 
field, A{K) = Xm=i 2^a{K)^ A " 1 -In the absence of a magnetic field we have Xm=i 2 3^0, A(K)a^ = and the 
fluctuating part is given by ,Xm=i 2 3-^o,A(K)a^ A \ In the momentum representation the model given in eq._Bl the 
external gauge field and without the static magnetic field takes form : S = Sq + S ex t , So is the free electron action 
in the presence of the Rashba term and S ex t is the electromagnetic C/ em (l) x SU(2) sp i n action. 

So = J dtj j^{* + (K,t)(-ihdt-E F )*(K,t) 

- y + (K,t)[e(K- J2 MK)<r {A) )+ J2 ^ A {K)o- {A) ]^{K,t)} (B2) 

A=l,2 A=l,2,3 

In momentum space the external vector potential A(q,t) and the scalar potential Ao(q,t) generate trough 
the A-C term a vector SU(2) gauge potential A(q) — Xm=i 2 A.A(q)o~( A ^ and a static SU(2) gauge potential 
J2a=i 2 3 Ao.a(<z)°' ( ^ • The external action [2£j is given by; S ext = S ( f xt ra " 1 + j : 



ext 



,4=1,2,3 



(B3) 



A=l,2 

We observe that the action is U(l) x SU(2) gauge invariant in momentum space. As a result the coordinate is 
SU(2) transformed , R{K) — > r(ir). Using Noether's theorem |30j | we obtain that the Spin - current is determined 
by the derivative with respect to the SU(2) gauge potential A^^iq)- 

- fray - 1 / |£* + («y<^- fi *(*) ( B5) 



where are the Pauli matrices. 

In order to study the electromagnetic response for an electric field without an orbital magnetic effect q x A(q) = 
we take the vector potential to be zero A(q) = and represent the electric field by a scalar potential A^ (q) which 

satisfies the condition E^? 1 ^ (q, t) = iq^A^^ (<f, t). (An alternative choice is to replace the vector potential by a time 

dependent vector potential. As a result we have to use a time dependent momentum (t) = (0) + iE[ ext ^t this 
approach has been used in the appendix A. ) From eq.-B6 we obtain that the spin velocity is given by the symmetric 
product between the spin operator and the velocity operator. Using the SU(2) transformation which diagonalizes the 

Rashba Hamiltonian introduced in chapter 777 with the velocity operator Vj(K) = ( dr j^^ )\t=o gives us according 

to eq.£?6 the spin velocity operator , ,l\ A '{K) — h{Vi(K), a^} + where {,}+ stands for the symmetric product and 

= W{K)u^U{K) represents the transformed Pauli matrix with A = 1,2,3. The transformed spin current 
operator in the Heisenberg picture is given by; 
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where = U + {K)a^U(K) with the Heisenberg representation of the matrix elements, Vi^ B (K ,t), r^(K,t) 
and ' (K ,t). The spin currents are computed by taking the expectation value with respect the Many-Body state, 
\F.S. >= \Kl,K F > , 

4 A) (q,t) =< F.S.\J$(q,t)\F-S. > (B8) 

Using Noether's theorem |3fj l3lj . we obtain that the spin currents are only covariantly conserved. 

d t J {A \q,t) + iq- J {A) (q,t) = 2e ABC [{^B + A , B (-£ i))J^ C) (q,t) + (lu b + A B (-q)) ■ J (C) {q,t)} (B9) 

where cabc^s the antisymmetric tensor which takes the values, 1,0,-1. We apply this equation to the Rashba model. 
The spin current components in the limit q — > obey the following continuity equation for the three polarizations: 
dt Jq 1 ' + *3 ' J^{o) = 0> ®tJo + *9 " J^il) — 0. Only the A = 3 component violates the continuity equation 

d t 4 3) + iq- fi 3 \q) = 2e 3BC LJ B (K) ■ / C) (g) = 2k so ( J x (1) (q) + j| 2) (q)) 
^ 2fc so (~) < F.5.| | ^*+(X,t)[i{^(#),a( 1 )} + + i{^(X),a( 2 )} + ]* ff (X,t)|F.5.> 



(27T) 

2fc so (^) | |^ <^.|vl/+(i?)[i{^ 3) (i?,t),(^(i?,t) x 

+ ^{^(K.i),^^,*).^)}^^)!^. > (BIO) 
z |K| 

The first term on the right hand side of the equation vanishes as a result of to the momentum integration .The second 
term vanishes as a result of the expectation value with respect the Pauli matrices. The external electric field induces 
a first order contribution to the spin velocity which vanishes after the momentum integration is performed. Therefore 
we conclude that the continuity equation is satisfied effectively. 

APPENDIX C: THE SPIN AND CHARGE HALL EFFECT IN THE PRESENCE OF A ZEEMAN 

INTERACTION 

In this appendix we will consider the charge Hall effect in the static linear response approximation. We will consider 
the current for the following cases : No orbital magnetic contribution, magnetic orbital current without spin Orbit 
and at the end we compute the charge Hall in the presence of the spin Orbit interaction and a weak magnetic field. 

In the absence of the orbital motion the charge Hall conductivity is given by; 

|^ [fF.D.(e so + b 3 - E F ) - f F . D .(-b 3 + e so - E F )] (CI) 

This result follows directly from equation 5 once we replace the Pauli matrix with the identity matrix /. The 
conductivity is zero in the absence of the Zccman interaction and vanishes at zero temperature. 

-The orbital effects are investigated for a periodic potential W(q) with the Bloch cigenfunctions u nj K(q)- n is the 

band index and q is the coordinate in the unit cell. The coordinate R is a function of q, R = q + f= q + i-§^- 

The curvature is determined by the band connection A rlym (K) = if u* m K (q)J^u n _ K (q). In the one band 

approximation the commutator is given by, [r^\ r^~\ = iQ(K, n) with the band curvature, £l(K, n) = x A n ^ n (K ). 

We apply this formalism to electrons in a constant magnetic field b perpendicular to the two dimensional plane. 
(For a strong magnetic field which is such that the magnetic flux per unit cell is a rational multiple ^ of the the flux 
unit $0 we obtain a magnetic unit cell with I cites and a reduced Brillouin zone. The complex eigenfunction u n ^K{q) 
has zero's in the Brillouin zone which is the origin for the U{\) Berry phase connections.) In order to study orbital 

effects we replace the Rashba Hamiltonian equation by, ho(K,r) = j—(K — jr-b x r — k so (a x e 3 )) 2 . We define the 

kinetic crystal momentum k which replaces the momentum K by, k = K — -^—b x r. The kinetic crystal momentum 
components do not commute in a magnetic field, [ki, k&] = —i-S-b. 

In the absence of the spin — orbit interaction the Bloch function in a magnetic field |32j gives rise to a U(l) Berry 

phase connection A n _ n (K) —if u* n K(o)w! u n,K{Q) HH and a Berry curvature Q(K,n) — x A n , n (K). In 



15 



the one band approximation we find [r^,r^l = ifl(K,n). We introduce the 4x4 Symplectic matrix to describe 
the commutators of the Cartesian coordinates and kinetic momentum, Jjj = where £ = (ri, T2, K%, K2) = 

£2, £3, £4)- The Heisenberg equation of motion are given by, 

(in) gjjg) =E<7 QMKAK)) (C2) 

dt ^ l ' 3 d ^{K) K ' 

We obtain the same equation of motion as the one given in ref. 34| (see eqs.l, 2). Using the equation of motion we 
compute the charge Hall current in the i = 2 direction, 

4%au = ^J2Tr[^^F F . D XEn(K),E F )] (C3) 

K,n 

"Tr" stands for the trace over the spin degrees of freedom. Fp,D,{K,Ep) = S(a,1)f FD (E n (K), Ep) + 5(<r,l 
)f F D (E n (K), Ep) is the Fermi-Dirac function and E n (K) are the Bloch eigenvalues as a function of the magnetic 
band index n. We find that the conductivity is given by, 

2 /> 

T / d 2 Kj2^(K,n)[f F . D (E n (K) +b 3 - Ep) + f F . D (E n (K) - b 3 - E F )] (C4) 

This result is in agreement with the result in ref. ( see eq.3). For strong magnetic fields the commutator 
M 1 ), r^] = ifl(K) describes the vorticity for the Hall wave function [24|. We find that the conductance is determined 
by the integral over the curvature and is given by the first Chern number [24| (see eqs.3.9 and 4.8). 

Next we consider the combined effect of a spin Orbit interaction and a weak magnetic field. For each Bloch band 
with magnetic eigenfunctions u n ^K(<T) we compute the following SU(2) matrix elements, (e{n)) n ,n + (<? ■ &)n,n where 
($l)n,n = [{— h m* Q K 2, - ° K\, \ga B {B + &))]«,« is a fictitious magnetic field for the band n and {e{K)) n ,n = E n (n) 
is the eigenvalue of the Bloch state and b 3 = j9 u b{B + b) is the total Zeeman field. We introduce an SU(2) 
transformation to rotate the Pauli matrix a in the direction of the fictitious magnetic field (3?) n ,n where {K) n ,n = 

~ 2lc^ x (f)n,n = K — 2lc x A n , n (K). For weak magnetic field we find that the curvature for the band n is 
given by, [r' 1 ), n 2 '] = H • £l{K,ri) + | [—a 1 sm-d(nn.n) + & 3 cos i?(/5 n ,n)] 27r5 2 ((«;) n! „). The first term is the Berry 

curvature of the magnetic Bloch function (generated by the connection) A n>n (K) — i J ,^ff a u^ K {q) ^u n> K(q) [33| . 
The second part is the curvature matrix due to the ££7(2) spin-orbit interaction. This part is obtained by replacing 
K — > (/t)n,n = K — -S-b x (r) n .n = K — ^-bx A njl (K). As a result we find that the Hall conductivity has two parts, 
a band contribution and a spin Orbit part; 

77 / d 2 Kj2^(K,n)[fF. D (E n (K) + e so + b 3 -Ep) + fF. D (E n (K) + e so -b 3 -Ep)} 
e 2 v 

+ 7^ $}/f.dOE*(0) + e so + b 3 - Ep) - f F . D (E n (0) + e so - b 3 - E F )\ (C5) 



In the present case the spin orbit contribution to the current (the second term) is determined by the magnetic 
cyclotron eigenvalues E n (0) » i£z{ n + h)> n = 0, 1) 2, 3, .... The spin Orbit contribution to the current decreases with 
the increase of the magnetic field faster than in the case where orbital effects are absent. 
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